Conditional Limit Results for Type I Polar Distributions 



Enkelejd Hashorva 



OO 



Department of Mathematical Statistics and Actuarial Science 
University of Bern, Sidlerstrasse 5 
CH-3012 Bern, Switzerland 
October 8, 2008 

Abstract: Let (Si, S2) = (-Rcos(0), i?sin(0)) be a bivariate random vector with associated random radius R which 
has distribution function F being further independent of the random angle 0. In this paper we investigate the 



' asymptotic behaviour of the conditional survivor probability ^ p>u (y) := P{pS\ + y 1 — p 2 S2 > y\Si > u},p <G 

. (—1, M when u approaches the upper endpoint of F. On the density function of O we require a certain local 

' asymptotic behaviour at 0, whereas for F we require that it belongs to the Gumbel max-domain of attraction. The 

q ' main result of this contribution is an asymptotic expansion of ^ p , u , which is then utilised to construct two estimators 

. for the conditional distribution function 1 — ^ PlU - Further, we allow O to depend on u. 

qq \ Key words and phrases: Polar distributions; Elliptical distributions; Gumbel max-domain of attraction; conditional 
limit theorem; tail asymptotics; estimation of conditional distribution. 

h' 

V\' 1 Motivation 

I Let (SijS^) be a spherical bivariate random vector with associated random radius R > (almost surely) with 
distribution function F. The random vector (X, Y) with stochastic representation 



(X,Y) = (Si, P S 1 + ^l-p 2 S 2 ), pe(-l,l) 



is an elliptical random vector (= stands for equality of the distribution functions) . If F is in the Gumbel max-domain 



lim P{Z U , P >pu + yy/u/w(u)\ = p\ Z > y/ y/l - p 2 }, VyeR, (2) 



\ of attraction with positive scaling function w, i.e., 

1 - F(u + x/w(u)) , N w _ 

u\x F 1 — t (U ) 

_J 1 where xf € (0, 00] is the upper endpoint of F, then Theorem 4.1 in Berman (1983) implies the following Gaussian 
, approximation 

> 

• rH 

X 

■ with Z u ^ p = Y\X > u and Z a standard Gaussian random variable (mean and variance 1). 

Berman's result shows that the Gumbel max-domain of attraction assumption is crucial for the derivation of ([2]). 
Conditional limit results for F in the Wcibull max-domain of attraction and (X, Y) a bivariate elliptical random 
vector are obtained in Berman (1992), Hashorva (2007b). The case F is in the Frechet max-domain of attraction is 
simpler to deal with, see Berman (1992). 

As shown in Cambanis et al. (1981) we have the following stochastic representation 

(Si, 5 a ) = (i?cos(e),i?sin(0)), (3) 

with R independent of the random angle O which is uniformly distributed on (— 7r,7r), i.e., (cos(O)) 2 possesses the 
Beta distribution with parameters 1/2, 1/2. 

When (cos(6)) 2 is Beta distributed, then the random vector (Si,S2) is a generalised symmetrised Dirichlct random 
vector. Generalisation of ((2|) for such (SijSa) is presented in Hashorva (2008c) with limit random variable Z being 
Gamma distributed (see below Example 1). 
Three natural questions arise: 

a) What is the adequate approximation of the conditional survivor function P{Z U . P > y} if O e (— tt.tt) is some 
general random angle with unknown distribution function? 

b) What can be said about the limit random variable Z? 

c) Does Z has a more general distribution if the random angle O = Q u varies with u! 
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In this paper we show that if Q u possesses a positive density function h u with a certain local asymptotic behaviour 
at 0, then we can answer both questions raised above. The generalisation of @ for bivariate polar random vectors 
(see Definition [T] below) satisfying |T]) is given in Section 3. Two applications of our results are presented in Section 
4. The first one concerns the asymptotic behaviour of survivor function of bivariate polar random vectors. In the 
second application we discuss the estimation of the conditional distribution function P{Z UtP > y}. Proofs and related 
results are relegated to Section 5. 



2 Preliminaries 

We shall explain first the meaning of some notation, and then we introduce the class of bivariate polar random vectors. 
A set of assumptions needed to derive the main results of this paper concludes this section. 

If X is a random variable with distribution function H this will be alternatively denoted by X ~ H . When H 
possesses the density function h we write X n h. 

In the following ip is a positive measurable function such that for all z £ (0, oo) 

V>0) < Kmax(z x \z X2 ), K > 0, A, £ (-1/2, oo), i=l,2, (4) 
where E{ip(W 2 /2)} > with W ~ $. Since E{ip(W 2 /2)} < oo we can define a distribution function <I> on_ZR by 

Loo exp(-s 2 /2)?/>(s 2 /2) ds 

We denote by \t„^, a, f3 > the Gamma distribution with density function x Q_1 exp(— /3x)/3"/r(a), x € (0, oo), where 
r(-) is the Gamma function. 

Next, we introduce the class of bivariate polar random vectors. Throughout the paper R denotes a positive random 
radius with distribution function F independent of the random angle £ (— 7r, 7r), and (Si, S2) is a bivariate random 
vector with representation j3]). In the special case is uniformly distributed on (— 7r,7r) for any two constants 0,1,0,2 
(see Lemma 6.1 in Berman (1983)) we have 



cu Si + 02^2 = yaj + a 2 2 Sx = J a 2 + a 2 5 2 , (6) 

hence linear combinations of spherical random vectors (i.e. the elliptical random vectors) are very tractable. 
If the random angle is not uniformly distributed on (— 7r, 7t), then ([6|) does not hold in general. In this paper we do 
not make specific distributional assumptions on 0. We assume however that the random angle possesses a positive 
density function h on (— 7T, 7t). 

Definition 1. A bivariate random vector (X,Y) is referred to as a bivariate polar random vector with coefficients 
Oi,bi,i = 1,2 if it has the stochastic representation 

{X,Y)^(a 1 S 1 + aaS 2 ,b 1 Si + b 2 Sa), (Si, S 2 ) = (#cos(0), iisin(0)), (7) 
where R ~ F and R > (almost surely) being independent of the random angle £ (—ir,ir). 

Clearly, bivariate elliptical random vectors are included in the above class, which is defined in terms of three compo- 
nents, a) the distribution of the associated random radius R, b) the distribution function of the random angle 0, and 
c) the deterministic coefficients 01, 0,2, bi, 62- In this paper we consider for simplicity the case 



ai = 1, a 2 = 0, and 61 = p,b% = yl — p 2 , p £ (—1, 1). 

We refer to p as the pseudo-correlation coefficient, and call (X,Y) simply a bivariate polar random vector with 
pseudo-correlation coefficient p. We have thus the stochastic representation 

(X, Y) = (Si, pSi + Vl - P 2 Sz), (Si,S 2 ) = (Rcos(Q),Rsm(e)), R~F, (8) 
with R > independent of 0. 

We note in passing that Si, S2 are in general dependent random variables. If Si and S2 are independent, for instance 
if R 2 is chi-squarcd distributed with 2 degrees of freedom and is uniformly distributed on (— ir,ir), then (X,Y) is 
a linear combination of independent Gaussian random variables. 
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Next, we formulate three assumptions needed in this paper: 
Al. [Gumbel max-domain of attraction] 

The distribution function F with upper endpoint xp is in the Gumbel max-domain of attraction satisfying |T]) with 
the scaling function w. Further, suppose that F(0) ~ and xf G (0, oo]. 

We formulate next an assumption for the second order approximation in JT]) initially suggested in Abdous et al. 
(2008). 

A2. [Second order approximation of F] 

Let F be a distribution function on [0, oo) satisfying Assumption Al. Suppose that there exist positive functions A, B 
such that 



1 - F(u + x/w(u)) 



1 - F(u) 



exp(-a;) < A(u)B(x) 



(9) 



holds for all u < xf large enough and any x G [0, oo). Furthermore we assume lim ll f. rF A(u) = 0, and B is locally 

bounded on finite intervals of [0, oo). 

A3. [Local approximation of h n ,n > 1 along t n ] 

Let h n : (—tt,tt) — > [0,oo),n > 1 be a sequence of density functions such that h n {6) = h n (~9),V9 G [0, 7r/2), and 
let t n ,n > 1 be positive constants tending to oo as n — > oo. Assume that for any sequence of positive measurable 
functions r n (s) = 1 + 0(s/t n ),n > 1, s > for all large n we have 



(10) 



h n (^Tn{s)\j^-\ = h n (l/y/t^)ip Tn (zT n (s)), Vs,zG[0,oo), 
where %j) Tn , n > 1 are positive measurable functions such that 

Vv„(s) — > ^(s), n —> 00 

locally uniformly for s,z G [0, oo) with -0 Tn satisfying ^ for all large n and all s G [0, et„) with e a fixed positive 
constant. 

Next we impose an assumption on the second order asymptotic behaviour of h n , n > 1 at 0. 
A4. [Second order approximation of h n , n > 1 along t n ] 

Suppose that Assumption A3 holds for some given sequence t n ,n > 1, and further for any sequence of functions 
T n (s) = 1 + 0(s/t n ), n > 1, s > for all large n we have 
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h n {l/\/t7i) 



ip(z) <a(t n )b n (z), Vs,zG[0,oo) 



(11) 



where a,b n ,n > 1 are positive measurable functions such that 

lim a(t n ) = 0, lim b n (s) = b(s), 

n — ^oo n — 'oo 

and b n ,n > 1 satisfy ^ for all n large. 



3 Main Results 

In this section we consider a bivariate polar random vector (X, Y) with pseudo-correlation p and representation (jSJ) . 
We are interested in the asymptotic behaviour of the conditional distribution Y\X > u n when u n tends (n — > oo) 
to the upper endpoint xp of F. Several authors have dealt with such conditional probabilities and their statistical 
estimation, see e.g., Gale (1980), Eddy and Gale (1981), Berman (1982, 1983, 1992), Heffernan and Tawn (2004), 
Abdous ct al. (2005), Heffernan and Rcsnick (2007), Abdous ct al. (2008) and Hashorva (2008b,c). Statistical 
modelling of conditional distributions is treated in the excellent monograph Thomas and Reiss (2007). 
The main assumption imposed on F is that it satisfies Assumption Al with the scaling function w. Such polar random 
vectors are referred to alternatively as Type I polar random vectors. The scaling function w possesses two crucial 
asymptotic properties: a) uniformly on the compact sets of J? 

lim ™(" + *M")) = 1, (12) 

U^Xp W(u) 
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and b) 



lim uw(u) = oo, lira w{u)(xp — u) = oo if xp < oo. (13) 



Refer to Falk et al. (2004) or Resnick (2008) for more details on the Gumbel max-domain of attraction. 

We derive in the next theorem the asymptotic behaviour of -Rcos(0„), with 0„ a random angle depending on n. 
Theorem 1. Let R be a positive random radius with distribution function F independent of the random angle 
n h n ,n > 1. Let u n ,n > 1 be constants such that u n < xf,u > 1 and lim^—Hx, u n — xp with xp G (0, oo] the 
upper endpoint of F. If F satisfies Assumption Al, and the density functions h n ,n > 1 satisfy Assumption A3 along 
t n := u n w(u n ), n > 1 with ip, ip Tn , n > 1, then we have 

poo 

P{R cos(0„) > u n } = {l + o(l))t-^ 2 h n (l/Vui)[l-F{u n )} exp{-x 2 /2)^{x 2 /2)dx, n -> oo. (14) 



If 0„ = 0,Vn > 1 not depending on n, then i?cos(0) has distribution function in the Gumbel max-domain of 
attraction with the scaling function w. Furthermore, the convergence in probability 

q n \Rcos(Q) — u n \ Rcos(Q) > u n A 0, n — > oo (15) 

holds for any sequence q n ,n > 1 such that linij^oo w(u n )/q n = oo. 

We note in passing that (jT4j) is obtained in Theorem 12.3.1 of Berman (1992) assuming that (cos(0„)) 2 is Beta 
distributed with positive parameters a, b. See also Tang (2006, 2008) for some important results on tail asymptotics 
of products of independent random variables. 
We state now the main result of this section. 

Theorem 2. Let (X n , Y n ), n > 1 be a bivariate polar random vector with representation |8"j). where p G (— 1, 1), R ~ F 
and 0„ i2 h n ,n > 1. Let u n ,n > 1 be a positive sequence such that u n < xp,n > 1 and limn^oo u n = xp. Suppose 
that F satisfies Assumption Al and h n ,n > 1 satisfy Assumption A3 along t n := u n w(u n ),n > 1 with ijj,ip Tn ,n > 1. 
If further limsup,^^ h n ((l + o(l)) / 'y/t^) j 'h n (l / \/t^) < oo, then for any x > 0, y £M we have 

lirn^ p|t„ < pu n + yu n /y/t^,X n < u n + x/w(u n ) X n > u„| = pjz < - p 2 , W < x|, (16) 

im£/i Z ^ ^ being independent of W ~ ^1,1, where ^ is defined in ([5]). 

Assumption A3 is somewhat cumbersome. If we consider random angles n not depending on n for all large n, a 
tractable condition on the local asymptotic behaviour of the density of 9„ is imposed below. 

Theorem 3. Under the setup of Theorem^ if Q n = Q i2 h, n > 1 and instead of Assumption A3 we suppose that the 
density function h of is regularly varying at with index 26 € ( — l,oo), then for any sequence u n < Xp,n > 1 such 
that linin^oo u n ~ xf we have 

P{X n >u n } = (1 + (l)) r( ^ 1 1 / ' 2) ^ 1/2 fe(yT7^)[l - F(u n )}, n^oo, (17) 

and X\ has distribution function in the max-domain of attraction of the Gumbel distribution with the scaling function 
w. Furthermore (| 1 5 j) holds for any sequence q n ,n > 1 such that linin^oo w(u n ) / q n — 00, and for x > 0, y G IR given 
constants (|16[) is satisfied with Z 2 ~ ^5+1/2, 1/27 an d % symmetric about independent ofW^ ^x,i- 
Wc present next an illustrating example. 

Example 1. [Kotz Type III Polar Random Vector] Let R ~ F be a random radius with tail asymptotic behaviour 

1-F(u) = (l + o(l))Ku N cxp(-ru s ), K > 0, 6 > 0, N G M, u^oo. (18) 

If ±2 ft, is a random angle independent of R we call (X, Y) with stochastic representation ([8]) a Kotz Type III polar 
random vector with pseudo-correlation p G (— 1, 1). If we set w(it) := r^it* 5-1 , it > 0, then 

,. P{R > u + x/w(u)} , 

lim — - — — r- 4 = exp(-x), WxGM 

u^oo P{R>U} 
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implying that F is in the Gumbel max-domain of attraction with the scaling function w. Suppose that h(9) = 
h(-0),V0 £ [0,tt/2), and further 

h{6) =c a , fc |sin(0)| 2Q - 1 |cos(0)| 2fc - 1 , 6e(-e,e), e £ (0,tt), 

where a, 6, c a .h are positive constants. Note that when 

1 T(a + 6) 
e-7r, c «,fc- 2r(a)r(6)' 

then (X, Y) is a generalised symmetrised Dirichlet random vector (see Hashorva (2008c)). It follows that Assumption 
A3 is satisfied with 

Ml/VQ = (l + o(l))Ca,6*i /3 " , V>00 = (2s) a - 1/2 , s>0, i„^oo 
and ft. is regularly varying at with index 2a — 1. By (jT4j) for u„ — > oo we have 

P{X > u n } = (1 + o(l))c a , 6 JC(2/(r<5)) a r(a)^- a5 exp(-r«*). 
Next, for any x > 0,y Theorem [2] implies 

lim p\y < pu n + yu]T Sf2 ,X < u n + xulr 5 X > u n \ = P{Z < yy/r5/{l - p 2 ), W < rSx}, 

with Z symmetric about independent of W ~ and Z 2 ~ 4 , x/2- Remark that if a = 1/2, then Z is a standard 
Gaussian random vector. When also b = 1/2, then (X, Y) is an elliptical random vector with pseudo-correlation p. 
In the next theorem we show a second order correction for the conditional limit result obtained in (|16[) which is of 
some interest for statistical applications. 

Theorem 4. Under the assumptions and the notation of Theorem^ if furthermore Assumptions A2 and A4 are sat- 
isfied where xp = oo and p £ [0, 1) , then we have locally uniformly for any z £JR (set z n _ p := pu n + zu n yj\ — p 2 /y/tn) 



P{y n > z n .p X n > u n \ = l-$(z) + ——JL=y'(z)+o(A(u n )+a{t n ) + —) ) n -> oo, (19) 
' \Jt n y/\ — p z \ t n y 

provided that 



max 



^ B(s)ds, J B{s) max(s Al , s X2 )ds^j < oo, 



where A, £ (— l/2,oo),i = 1,2 are i/ie constants related to Assumption A3. 

Remark 1. a) Abdous et al. (2008) show several examples of distribution functions F satisfying Assumption A2. 
The assumptions on h can be easily checked for common distribution functions using Taylor expansion, 
b) If we assume h is regularly varying with index 25 £ (— 1, oo) instead of the Assumption A3 and modifying A4 
accordingly, then (fT9"j) holds with ^ := ^s+1/2.1/2: provided that 



max 



J B(s)ds,j B(s)s s ds j < oo. 



4 Applications 



In this section we present two applications of our asymptotic results: a) we obtain an asymptotic expansion for 
the joint survivor probability of polar random vectors, and b) we discuss briefly the estimation of the conditional 
distributions of such vectors. 



4.1 Tail Asymptotics 

Let (X,Y) be a bivariate polar random vector with pseudo-correlation coefficient p £ (—1,1). Assume that the 
distribution function F of the random radius R has an infinite upper endpoint. In various situations quantification 
of the asymptotics of the joint survivor probability P{X > x, Y > y} is of interest when x, y become large. Our 
asymptotic results in Section 3 imply an asymptotic expansion of this survivor probability, provided that (X, Y) 
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is of Type I. Explicitly, under the assumptions of Theorem [3] we obtain for any x > 0, y G IR and u large (set 
x u := u + x/w(u),y u ^ p := pu + y^u/w(u), u > 0) 

P{X>x u ,Y >y u , p } = {l + o(l))exp(-x)[l-y s+1/2A/2 (y)]P{X >u}, u -> oo. 

In our asymptotic result the sequence y u ^ p increases like pu since by (13) 

y 

y u , P = (H , )pu = (1 + o(l))pu, u^oo. 

p\Juw(u) 

It is of some interest to consider also constants y UtP — cu,u > 0, with c 6 (p, 1]. In view of Theorem [3] for any 
c G (— oo,p) we have 

P{A > a; u ,Y~ > y u , c } = (1 + o(l)) exp(-a;)P{X > u}, u -> oo. 
When c G (p, 1] the joint survivor probability P{X > x u , Y > y u . c } diminishes faster than P{A > it}, i.e., 

Um P{X>x u ,Y>y u , c} = Q 

n->oo P{A > u} 

If (X,Y) is a bivariate elliptical random vector we may write (see Hashorva (2007c)) 

P{X>u,Y>cu} = (l + (l)) ap '^ Ap ' cl ~ f(a ^ ) , u^oo (20) 

27r uw(a p>c u) 

for any c G (p, 1] with 

a p , c := v/(l - 2cp + p 2 )/(l - p 2 ) G (1, oo), tf„, c := t^^A £ (0,oo). 
In a forthcoming paper we extend (|20[) to the case of Type I bivariate polar random vectors. 
4.2 Estimation of Conditional Distributions 

Let (Xi,Yi),i < n,n > 1 be independent and identically distributed bivariate polar random vectors with pseudo- 
correlation coefficient p G (— 1, 1) and random radius R ^ F. Define the conditional distribution function 

* p , x {y) := P{Yi < y\X x > x}, x,y eM. 

Suppose that xf — oo and F satisfies Assumption Al. As in the elliptical setup (Abdous et al. (2008)) also in the 
general case of polar random vectors estimation of the conditional distribution function & PtX can be motivated by our 
novel limit results, since under the assumptions of Thcorcm[3] we have (set t u := uw(u), u > 0) 



sup 

yelR. 



^ P A u [p + v\fiTQ) - *5+i/2,i/ 2 (y/v / i - p 2 ) ->o, u^oo, (21) 



with 28 the index of the regular variation of h at 0. Under Assumptions A2 and A4 we obtain additionally the second 
order asymptotic expansion 

$p, u (u[p + l/t u + yy/TJQ) = * 4 +i/2,i/2(j//Vl - P 2 ) + 0{A{u) + a(u) + i), u -> oo. (22) 

These approximations motivate the following estimators of 3> p x for x large and y positive, namely 



(v) v-P" x ) 



n > 1, 



and 



, \ t /"V — Pn(a; + l/w)fi(x))\ 

* S """ !=W ^( V(i-«)»M.w )' n>1 ' 



where p„ is an estimator of p, and is an estimator of w(-). 



7 



An estimator of p n can be constructed considering the relation between p and the expectation E{Y}, provided that 
the latter exists. Estimation of 6 and w are difficult tasks. If the scaling function w (related to the Gumbel max- 
domain of attraction of F) is simple, say w(u) — cyu 1 , c, 7 > 0, u > 0, then an estimator w n is constructed by 
estimating separately c and 7 from X\ , . . . , X n (recall X\ has distribution function in the Gumbel max-domain of 
attraction with the scaling function w). See Abdous et al. (2008), Hashorva (2008a) for more details. 
In practical situations also the constant 5 might be unknown and therefore has to be estimated. One possibility of 
estimating S is to utilise fTTJl . 

We note that for elliptical random vectors <5 = 2 and both estimators ^ p ]l t n and ^f p %, n are suggested in Abdous et 
al. (2008). Since we estimate both c and 7 from X\, . . . , X n and not from the observations of the random radius R, 
our estimators above differ from those in the aforementioned paper. 

5 Related Results and Proofs 

Set in the following 

a P (x,y) := s/l + {{y/x) - p) 2 /(l - p 2 ), a*(x, y) := a p (x,y)x/y, x,yeM,y^0, pE (-1,1). (23) 
For 1 < a < b < 00, x > constants, and h, F two positive measurable functions we define 

f b 1 

J(a,b,x,h):= / [1 - F(xt)]h(t) dt. (24) 

J a t\/t 2 - 1 

If b = 00 write simply J(a,x,h) suppressing the second argument. Write h(-) and h p (-) instead of /i(arccos(l/-) and 
/i(arcsin(l/-) — arcsin(p)), respectively. 

Next, we shall prove two lemmas, and then proceed with the proof of the main results. The first lemma is formulated 
for F with infinite upper endpoint. It generalises Lemma 5 in Hashorva (2008b) for bivariate elliptical random vectors. 
If F has a finite upper endpoint, say xp = 1, then a similar result holds. Statement b) and c) should be reformulated 
requiring additionally that x 2 + 2pxy + y 2 < 1 — p 2 with |y| £ [0,1]. 

Lemma 5. Let the random radius R <~ F be independent of the random angle S (— n,ir) and define a bivariate 
polar random vector (X,Y) with pseudo- correlation p G (— 1, 1) via (|8|). If the upper endpoint xf of F is infinite and 
Q possesses a density function h such that h(9) = h(—8),6 £ [0, 7r/2), then we have: 

a) For any x > 

P{X>x} = 2J(l,x,h). (25) 

b) For any x > 0, y G (0, x] such that y/x > p 

P{X>x,Y>y} = J(a p (x,y),x,h) + J(a*(x,y),y,h p ). (26) 

c) For any x > and y/x G (0, p), p > 

P{X>x,Y>y} = 2J(l,x,h)-J(a p (x,y),x,h) + J(a* p (x,y),y,h p ). (27) 

Proof. Since the associated random radius R is almost surely positive being further independent of and h(—ff) = 
h(6), cos(— 6) = co&(9),M9 G [0, ir/2) for any x > we obtain 

,tt/2 

P{X > x} = 2 P{R > x/cos(9)}h(6)d9 
Jo 

= 2 Tll-Fixs)} ^fL= ds = 2J(l,x,h). 
J 1 sv s — 1 

We prove next the second statement. By the assumptions (X,Y) = (i?cos(0), i?sin(0 + arcsin(p))). Consequently 
for x > 0, y > two positive constants 

P{Si > x, pSi + ^Jl - p 2 S 2 > y] = P{R cos(0) > x, i?sin(0 + arcsin(p)) > y}. 
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Since sin(arcsin(p) + &)/ cos((9) is strictly increasing in 9 G [— arcsin(p), ir/2] with inverse arctan((- — p)J -J 1 — p 2 ) (see 
Kliippelberg et al. (2007)) we have 

,7T/2 

P{X>x,Y>y} = / P{R > x/ cos(6)}dQ(6) 

J arctan((i//a;-p)/ y/l-p 2 ) 



P{R > yj sin(0 + arcsin(p))} dQ((9), 

- arcsinfo) 



r arctan((y/a;-p)/ A /l-p 2 ) 
- arcsin(p) 

with Q the distribution function of 0. Transforming the variables we obtain for y/x > p 

P{X>x,Y>y} = J(a p (x, y), x, h) + J(a*(x,y),y,h p ), 
and if y/x < p with x, y positive 

P{X>x,Y>y} = 2J(l,x,h) - J(a p (x,y),x,h) + J(a*(x,y),y,h p ), 

hence the proof is complete. □ 

Lemma 6. Let F be a distribution function with upper endpoint ij? £ (0,oo] satisfying further ([I]) with the scaling 
function w and let 1 < a n < b n ,~f n > l,u n £ (0, x.p),i n := u n w{^ n u n ),n > 1 be positive constants such that 

b n u n < x F ,n > 1, lim 7„ = 7 G [l,oo), lim 7„u„ = lim b n u n = x F , (28) 

n — >oc n — >oo n — >oo 

and further 

lim i n (a„ - 7„) = £ G [0, 00), lim t n (b n - j n ) = 77 e [£, 00]. (29) 

n — >oc n — >oo 

Lei h,r,Tp n ,n > 1 be positive measurable functions. Assume that for some e > 

Hln + s/t n ) = r(7„,i n )V„(s), VsG[0,rf„] (30) 

and 

"0n(s) — * VK S ) € [0,Oo), 71 — > CO (31) 

locally uniformly with ip n satisfying (|4]) /or all n > 1, s G [0, ei n ] wii/i Aj,i = 1,2 £ (c, 00). Suppose further 
/~ fe(s)(sVs 2 - 1) _1 ds < A' < 00, Vn > 1. 
7f 7 G (1, 00) and c = — 1 

J(a n ,b n ,u n , h) = (l + o(l)) — -= = — — — f exp(— x)ip(x) dx, n—>oo. (32) 

7V7 2 — 1 tn j£ 

b) When 7 = 1 and lim^oo t n (j n — 1) = r G [0, 00), i/ien 

J(a n ,b n ,u n ,h) = (1 + o(l))r(7„,t„) — - / exp(-x) ^>(x) dec n — ► 00, (33) 



/*„ J ( y/2x + 2t 

provided that c = — 1/2 i/£ = T = and c = — 1 otherwise. 
Proof. Set in the following for n > 1 

< : = Tn^n, i« := u n w«), ^„(x) = 7„ + x/i„, V„W : = 77-^i==iW. a; > 

ww^w - 1 

and £„ := t n (a„ - 7„), r?„ := t„ (b n - 7„). Since lim rwoo u* = x F , then ((13J) implies lim„_ >o0 t n = oo, w(u* )(ijr - u* ) 
co. If ci, C2 are two arbitrary constants such that C2 > c\ > 1 for all n large we have 

1 Z 100 1 

[l-FM]/i(s)^==& < [1-AK( 7 + C2 ))] / h(a)—==da 
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and 



f n 1 
J(a n7 b n ,u n ,h) > / [1 - F(u n s)]h(s) — ds 

Ja n S\/S- - 1 

> f \l-F{u n s))h{s) 1—r ds 

> [1 - F(«n(7 + ci))] / 7+ 1 ft(g) r^- ds. 

J a n S\ S 1 

Assume that xp = oo. Since 1 — F is rapidly varying (see e.g., Resnick (2008)) i.e., 

1 - F(u n x) 
llm ~i Wi — r = °' Vx > 1 



for any e* > we obtain 

p+e" ^ 

J(a n ,b n ,u n ,h) = (l + o(l))/ [1 - F(«„s)]/i(s) — ds. 

Ja„ svs - 1 

If 7 G (1, oo), then 

locally uniformly for s > and ^* satisfying ([4]) for all s <G [0, £t n ),s > 0. As in the proof of Lemma 7 of Hashorva 
(2007a) for any e > as n — > oo wc obtain 

J(a n ,b n ,u n ,h) 

f .min(7„+T; 7l /t 7l ,7+e) ^ 



(1 + o(l)) / [1 - F(u„ a )]ft(i)— ===== d* 

fl+nfl^ / ,min (')«.*«('r-7>.+e)) 1 

[1 - F« + x / w (u* n )))h(l n (x))—-^=== dx 



tn j£ n l n (x)y/l*(x) - 1 

(1 + o(l)) r(7„, t n ) /■ min ('?"' t »(T-T™+ e )) 

[1 - F(u n + ar/«;«))]^(a;) da; 
_L=%M[i_ f(<)] J\xp(~xMx)dx. 



7v 7 2 — 1 ^™ 
( 1 

7 V 7 2 — 1 ^« 



Next, if 7 = 1 redefine 



<(s) := -= —. —, ===Vn(s), n>l, S >0. 

V*n'n(s)V^(s) - 1 



We have 

<(«)-> / (S) =:V>*(s) 
" V2r + 2s V ; 

locally uniformly for s > 0. Hence as in the proof above for e > and n-tcowe obtain 



1 



min(?7 n .et n ) 



J(a n ,b n ,u n ,h) = (i + o(l))-=/ [l-F(<+x/w(<))]V£(x)da: 



,'t n j£ 

1 - Ffit ) f v 1 
(l + o(l))r( 7 „,i„) -=L^ / exp(- a; )-===V'(x)dx. 



Similarly, the asymptotic results follow when x^? G (0, oo), hence the proof is complete. □ 
Proof of Theorem Q] We consider for simplicity only the case xf = oo. For all n large (|25|) implies 

/•oo 1 x 

P{R cos(0„) > u n } — 2 [1 — F(u„s)]/i ra (arccos(l/x)) dx. 

7i a: Vx 2 - 1 

We have (set t n := u n w(u n ), n > 1) 



arccos(l/(l + s/**)) = "T^ 1 + °( S A«)) = : \/2s/t n T n (s), n -> oo 

V tn 
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locally uniformly for s > 0. Hence the Assumption A3 on h n implies 



h n (arccos(l/(l + s/t n ))) = h n (T n {s)\j2s/t n ) = h n (l/y/t^)ip Tn (sT n (s)), s > 0. 

Applying Lemma [5] with r — 1 = 7 = 7 n = o n = 1, n > 1 and fe n = 00, n > 1 we obtain 

1 _ p( u ) r°° 

P{R cos(6„) > u n } = (l + o(l))MWQ 7=^ exp(-s 2 /2)^(s 2 /2)ds. 

\t n J —00 

If h n = h,n > 1, then by the Assumption A3 we have linin^oo h{l / y/tH) /h{y n / \ft^) = 1 for any sequence y n ,n > 1 
such that lim n ^oo y n = 1. Consequently the self- neglecting property of w in (|12|) implies 



P{i? cos(0) > u n + x/w(u n )} 1 - F(u n + x/w(u n )) 
hm — — = hm — r = exp(— x), vx eJR. 

n->oo P{R COS(B) > U n \ ra-»oo 1 — F{U n ) 

Hence for any z > 

r „ . , . „ , P\R cos(0) > u n + z/q„} 

P{q n \RcoB{e) - u n > z Rcos 6 > u n } = DfD — 

P{Rcos(Q) > Un + (z/w(u n ))(w(u n )/q n )} 
P{Rcos(Q) > u n } 

— > 0, n — > oo, 

thus the result follows. □ 
Proof of Theorem Set for n > 1 and 2 e.ff? 



fn = z^u n /w(u n ), Xn ■= v n /u n , a n := y/1 + (v n /u n ) 2 , t n := u n w(u n ), n > 1 

and write in the sequel h n {-) and h n , p (-) instead of /i„(arccos(l/-) and ft,„(arcsin(l/-) — arcsin(p)), respectively. 
Since linin^oo t n — 00 by the assumptions on h making use of Q13p and (|15p we retrieve the convergence in probability 

a/ w(u n )/u n (X n ~ u n ) X n > u n A 0, n — > 00. 

Consequently, it suffices to show the proof for p = 0. Next, we prove the convergence in distribution 

y/w(u n )/u n Y* — > Z ~ 71 — > 00, 

with y,* = T n |A„ > u n and $ defined in ([5]). Since = v n /u n > p = holds for all large n, we have in view of 
Lemma [3] for all large n 

P{X n > u n ,Y n > v n } = J (a n , u n , h) + J (Xn 
where a n = 1 + (1 + o(l))z 2 /(2i n ), n — > oo. As in the proof of Theorem [T] we obtain for the first term 
J(a n ,u n ,h) = {1 + °0))Kp-/y/Q {1 _ F{Un)] r exp{ - x y 2 )^(x 2 /2) dx. 

V <"n J v. 



Further, for any s > (set l n (s) := x n X + s/(v n w(u n )) we have 

TTT = — TT = V 2 /*«7T — 7T> 71 -> °°- 
Consequently, the assumption on h implies for all s > 

h n ,p(ln(s)) = h n (—j==T n (s))) = h n (l/y / t n ~)lp Tn (T n (s)z 2 /2) , 



where r n (s) := 1 + 0(s/t n ), s > 0, n > 1. Hence 

1 



in(«) V(^( fi )) 2 - 1 



VpGn(s)) = /i„(l/V^)x„ /2 ^r„(r„(s)z 2 /2), n 
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As in the proof of Lemma [6] we have thus 



JiXn 1 a n,Vn, h p ) 



[1-F(v n t)]h n , p (t) J_ 

Xn'Ct„ - 1 



dt 



h n (l/VQ 



1 - F(u n 
v„w(u n ) 



3/2 



1 - F(v n l n (s)) 

t„[a„-l] l--F'(Wn) 



(i + (i))/ i „(i/V^) L -^x^ / V(^ 2 /2) 

o(J(l,u„,/i„)), n — > oo 



VV„(t„(s)z 2 /2J ds 

1 - F(it„ + s/w(u n )) 



: /2(l+o(l)) 1 - ^(u„) 



implying 



lim P{Y* > z\J u n /w(u n )} 



lim 

n — >oo 



P{Y n > Zy/u n /w(u n ), X n > U n } 



P{X > U n } 

f™exp(-x 2 /2)ip(x 2 /2)dx 
/^ oo exp(-x 2 /2)V(x 2 /2)rfx 



i - #0). 



Thus the proof is complete. □ 
Proof of Theorem [3] By the assumption on h we have h(s) = s 2S L(s) for all s > in a neighbourhood of with 
L(s) a positive slowly varying function such that lim t _»o L(ts) / L(t) = l,Vs > 0. Furthermore, by Proposition B.1.10 
in de Haan and Ferreira (2006) we have for any e > 0, £ > 



h{t S ) 



h(t) 



< £ max( S 25 -«, S 25 +«) 



(34) 



holds for any s G (0,to(e,£)/t),t £ (0, 1) with to(e,£) some positive constant. Since for positive constants t n ,n> 1 
such that linin^oo t n = oo 



h(y/2 S /t n ) 



{2s) s =: ip(s), Vs > 



the result follows along the lines of the proof of Theorem [2] utilising further J34 
Proof of Theorem [4] Set forn > l,z £M and \p\ < 1 



□ 



v n := pu n + zy/i - p 2 \Ju n /w{u n ) 1 \n ■= V n /u n , t n := u n w{u n ), a„ := + (1/x 2 - p)/{l - p 2 )- 
In view of Lemma[5]for all large n we have P{X n > u n } = 2J(1, u n , h n ) and further for p > 0, z > 

1 



P{Y n > v n \X n > u n } 



2J(l,U n ,h n ) . 



J{oL n , u n , h n ) + J ( ^ n oz n: v n: h n p) 



In order to show the proof we need to approximate J(l + ^j—,u n , h n ) and J(x„ 1 Q!n, f n , h ntP ). We have 



a„ = 1 + — + 0(l/i 2 ), n>l. 



As in the proof of Theorem [5] we obtain 



J(a n ,u n ,h) = (1 + o(l))h n (l/y/tn) 



1 - F(Un) | 



1 - *(z)] / exp{-x 2 /2)tp(x 2 /2) dx, 



Assumptions A3 and A4 imply ip(x) < max(x Al , x x ' 2 ) and b(x) < m&x(x Xl ,x x ' 2 ) for some Aj, A* £ (— 1/2, oo). Conse- 
quently, the assumptions on B and b yield 

f°° 1 f 00 1 
exp(— x)ip(x)\/xdx < oo, / B(x)tp(x) . da: < oo, / cxp(— x)b(x) : dx < oo. 



V2 



2.7- 



Define for all x > and n > 1 



9n{x) 



(1 + x/t n )y/t n (l+ X /t n ) 2 -t n 
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For all x > 0, n > 1 wc have g n (x) < and further 



2.i- 



9n{x) 



< 



1 - 



1 



l + x/t n l + x/t r 



S2X y / 2x + X 2 /t r , 



< 5/4Vi/(t„V2). (35) 



Consequently, for £ > and n large we obtain 

A/ trt 



ft„(i/VQ[i-F((i„)] 

1 - F(«nt) 



J(l + Q/t n ,u n ,h n ) - / exp(—x)ip(x)—, = dx 

vzx 



-h n (t)- 



1+f/tn l - F(w„) /i„(l/VQ tx/F - T 7 C 
1 - F(u„ + x/w(u n )) h n (l + x/tn) 



dt — exp(—x)'il)(x)—;=dx 



< 



i - ^K) 

1 - F(u n + x/w(u n )) 



+ / exp(— x) 



1-F(u n ) 

h(l + x/t n ) 



h n (l/y/U) 
- exp(— x) 



1 



g n (x) dx — / exp(— .t)?/;(x)— p= dx 



2x 



h n (l + x/t n ) 



ip(x) 



fc»(i/VQ 

c/.c 



< A(u n ) / B(x)ip n (x)— = dx + I exp(-x)g n (x) 



dx 



h(l + x/t„) 



+ / exp(—x)ip(x) g n (x) 



1 



< A{u n ) I B{x)tp(x) dx + a(t n ) / exp(— x) b n (x) dx + 

V2x 



1 



2x 



i/j(x) 



5V2 
8t n 



dx 



exp(—x)ip(x)y / x dx 



Since 



we have 



= 0(.A{u n ) + a(t n ) + l/t„) =: R n (u n ). 



exp(—x)ip(x)—= dx 
z 2 /2+o(i/t„) v2x 



cxp(— x)ip(x)—= dx = 0(l/t n ), n — > oo 
z 2 /2 v 2 x 



ft„(l/VQ[l-%)] 



1 

exp(— x)ip(x)—== dx = R n (u n ). 
72 v2a; 



Assume for simplicity in the following that p > and z > 0. The other case can be established as in the proof of 
Theorem[3J We obtain the first order asymptotic expansion (set l n (s) := Xn 1 + s/(v n w(u n )), s > 0,rt > 1) 



[1 - F(v n t)]h n , p (t) J_ 
,u„/t-„ IVI - 1 



1 - FK) 
w n w(w„) 



„2 



1 - F(u„ + s/w(u n )) L 



V n w{u n ) J tn [a n -1] 1-F(u n ) 



^ TB («72(l + o(l)))ds 



1 - F(u„ + s/w(u n )) 



V n w(u n ) ^2/2+0(1/(2) 1 - F(u n ) 



i> Tn (zy2(l + o(l)))ds 



(l + o(l))-^=V(^/2)ft»(l/VQ 1 f ( " n) f° ex P (-.s)ds, n^oo 

V 1 - &n J 2 2/2 



Define next for n > 1 and s > 



We have for all s > 0, n > 1 



Sn(s) 



g n (s)-p < (—= + s/t n )K, g n (s) < pK, p := p 2 /^l- p 2 , 
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with K > 1 a positive constant. Next, for any £ > we may write 

v n w(u n ) 



h n {l/JQ[l-F{u n )]" x " n ^ n " n ' , " n ' p 
1 - F{vJ n (s)) h n , p {l n (s)) 
c 1 - F(u„) Ml/vO 

1 - F(u n + s/w(u n )) h„ iP {l n {s)) 



< 



< 



J(ln((),V„,h n>p ) -pV(C) ex P(-C) 

/>oo 

(s) ds — /5 / V(C) ex P( — s ) <^ s 
g„(s)ds- J o^(C)exp(-s) ds 



1-F(u„) hnO-fy/Q' 



1 - F(u„ + s/w(u n )) 



exp(— s) 



1 - FK) 

h n) p(l n (s)) 



exp(— s) 



h n ,p(ln(s)) 



h n {l/\/t^) 



g n (s) ds 



9n(s) - pip(0 



ds 



^(l/v^) 

) / B(s)ds + p J exp(-s) 

VpGn(s)) 



h n ,p(ln(s)) 



h n {l/\/t7i) 



ds 



+ / exp(-s) g n (s) - p 



ds 



< ^(u„)KpV(C) / S(s) ds + pa(t n )b n (0 exp(-C) + ^(C) / exp(-s) 5 „(s) - jo 



ds 



= Rn(u„)- 

Hence for all n large since a n = 1 + z 2 / (2i„) + 0(1/^) we may write 



^(1/^(1-%)) 
Consequently, as n — > oo 



and 



J(X« Oi n ,v ni h niP ) - pip(z /2)exp(-z /2) = R n (u n ). 



J(a n ,u n ,h n ) = h n (l/\/tn) 



1 - 



[1 - *(z)] / cxp(-a; 2 /2)^(.T 2 /2) cfe + i? n (w n ) 



i/;(z 2 /2)cxp(-z 2 /2)+i?„K) 



implying 



p{y„ > v n \x n > u n } = i-^(z) 



exp(-z 2 /2)^(z 2 /2) 



tn S_ oo eM-x 2 /2)^{x' 2 l2)dx 



+ R n (u n ), n — ► oo. 



The proof for z < follows with similar calculations, hence the result. □ 
Acknowledgement: I would like to thank two referees of Extremes for careful reading, several suggestions and 
corrections received in April 2008. 
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